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u  abstract.  ye  are  concerned  with  the  following  reliability  problem:  A  system  has 
k  different  types  of  components.  Associated  with  each  component  is  a  numerical 
value.  Let  { aJ }  (J  ■  l,...,k)  denote  the  set  of  numerical  values  of  the  k 
components.  Let  R(a  ,...,a*)  denote  the  probability  that  the  system  will 
isfactorily  fi.e.  Rfa1 , . . .  ,a*)  is  the  reliability  of  the  syst 


perform  sat! 


system) 

R(aA, . . . .a*)  has  the  properties  of  a  joint  cumulative  distribution 


and  assume 
function. 

Now  suppose  a|  <_. .  .<^  are  n  components  ol  .ype  j  (j  -  l,...,k). 
Then  n  systems  can  be  assembled  from  these  components.  Let  N 

denote  the  number  of  systems  that  perform  satisfactorily.  N  is  a  random  variab} 
whose  distribution  will  depend  on  the  way  the  n  systems  are  assembled.  Of  all 
different  ways  in  which  the  n  systems  can  be  assembled,  the  paper  shows  that 
EN  is  maximized  if  these  n  systems  have  reliabiliry  R(aj , . . .  ,aj[)  (i  »  1,.. 

The  method  used  here  is  an  extension  of  a  well  known  resulc  of  Hardy,  Llttlewood 
and  Polya  on  suras  of  products.  Furthermore,  under  certain  conditions,  the  same 
assembly  that  maximizes  EN  minimizes  the  variance  of  N. 

Finally,  for  a  similar  problem  in  reliability,  it  is  shown  that  for  a 
series  systems  a  construction  can  be  found  that  not  only  maximizes  the  expected 
number  of  functioning  modules  but  also  possesses  the  stronger  property  of 
maximizing  the  probability  that  the  number  of  functioning  modules  is  at  least 
r,  for  each  0  <  r  <  n. 
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ON  OPTIMAL  ASSEMBLY  OF  SYSTEMS 


Cyrus  Derman,  Gerald  J.  Lieberman  &  Sheldon  M,  Ross 


[0]  Summary 

We  «re  concerned  wlth^the  following  reliability  problem:  A  system 
has  k  different  types  of  components.  Associated  with  each  component 

‘  -P  * 

is  a  numerical  value.  Let  {aJ}  (J  ■  l,...,k)  denote  the  set  of 

v,  k  ' 

numerical  values  of  the  k  components.  Let  R(a  , ...,a  )  denote  the 

1  k 

probability  that  the  system  will  perform  satisfactorily  (i.e.  R(a  ,  ...,a  ) 

1  k 

is  the  reliability  of  the  system)  and  assume  R(a  ,.  ,a  )  has  the 
properties  of  a  joint  cumulative  distribution  function. 

Now  suppose  aj  are  n  components  of  type  j  (j  *  l,...,k). 

Then  n  systems  can  be  assembled  from  these  components.  Let  N  denote 

the  number  of  systems  that  perform  satisfactorily.  N  is  a  random  variable 

whose  distribution  will  depend  on  the  way  the  n  systems  are  assembled. 

Of  all  different  ways  in  which  the  n  systems  can  be  assembled,  the  paper 

shows  that  EN  is  maximized  if  these  n  systems  have  reliability 
1  k 

R(a^ . a^)  (i  =  l,...,n).  The  method  used  here  is  an  extension  of  a 

well  known  result  of  Hardy,  Littlewood,  and  Polya  on  sums  of  products. 
Furthermore,  under  certain  conditions,  the  same  assembly  that  maximizes 
EN  minimizes  the  variance  of  N. 

Finally,  for  a  similar  problem  in  reliability,  it  is  shown  that 
for  a  series  systems  a  construction  can  be  found  that  not  only 


1 


maximizes  Che  expected  number  of  functioning  modules  but  also  possesses 
the  stronger  property  of  maximizing  the  probability  that  the  number  of 
functioning  modules  is  at  least  r,  for  each  0  <  r  <  n. 


[1]  An  Optimal  Assignment  Theorem 

A  well-known  result  appearing  in  Hardy,  Littlewood,  and  Polya  [1] 

asserts  if  a.  <  a„  <...<  a  and  if  b,  <  b„  <...<  b  ,  then 

1—  2—  —  n  1  —  2  —  —  n 

n  n 

I  a.b.  >  J  a.b . ,.s 

i-l  1  1  -  i.i  i  <KD 

where  ip  is  any  permutation  of  the  integers  1,  2,...,n. 

Also  in  [1]  is  the  generalization  that  for  any  k  2  and  if 
0  <_  a“  . . .  <_  a“  ,  j  -  1, . . .  ,k  then 


k 

n 

j-i 


> 


I 


i-1 


a 


Vu 


where  if/^(i)  ■  i  (i  -  l,...,n)  and  i|»2, . . .  are  any  k  permutations 
of  1,  2,...,n.  Proof  of  both  results  follows  from  establishing  the 
inequalities  for  the  special  case  of  n  ■  2  and  then  resorting  to  a 
standard  argument  regarding  pairwise  interchanges.  First  let  n  ■  2 
and  k  =  2.  Then 


a^b^  +  ~  a1^2  ”  a2^1  “  a^Cbj^^)  -  ^(bj-bj)  ■  (a^-a2)  (b^-b2)  2.  0 

since  a^  <_  a2  and  b^  <_  bj.  Now  assume  the  result  is  true  for  n  ■  2 


2 


and  k  -  2,  3 . K-l.  Then  consider 


i  ^  «j  (1)- 
i-i  j-i  V  ' 


If  for  each  J  ■  2,...,K,  (1)  -  2,  (2)  ■  1.  Then 


II  a“  ,,v-  b0 
J-2  *j(1)  2 


(say) 


II  a“  -  b, 

j-2  *j(2)  1 


1  b2  * 


(say) 


and,  hence. 


K  K 

n  a  +  n  aJ 
j-1  j-1 


l,  i  1, 

-  a1b1  +  a2b2 


i  ajb2  +  a2b1 


K  K 

n  a,  +  n  a, 
j-i+j(i)  hV8 


If,  for  at  least  one  j  (say,  j  -  2),  (1)  -  (2)  -  2,  then  by 

the  induction  assumption 


a?  fa?  If  J  J  +  fa?  a?  5  a? 


1  l  1  j-3  *j  U)|  T  L“2  2  ,:3  ‘♦j  <2) 


1^*1  11^1 
a.  n  a”  +  JI  a"  +  (a„-a,)  n  a- 

1  ij-l  V11  j-2  2  1  Lj-2  *j(2)J 


1^1  11  ^  j 

<.  a1  II  a^  +  II  a^  +  (a_-a.)  II  a** 

Lj-2  j-2  lj-2  Vj 


3 


±  ai 


K  K  ' 

n  aj  +  n  aj 

J-2  J-2 


(a,  - 


n  a, 
j-2 


n  a 

j-1 


J 


n  a: 
J-1 


j 


Hence  the  Inequality  holds  for  n  -  2,  and  all  integer  values  of  k. 

Returning  to  arbitrary  n,  if  1 1>2 . ^  are  not  all  equal  to 

then  there  exist  two  values  of  1  (say  1*1,  and  1-2)  such 

that  ’J'j(l)  L  $j(2)  does  not  hold  for  all  j  -  2 . k.  However,  by 

considering  the  permutations  (j  *  2 , . . .  ,k)  which  are  the  same  as 

'J'j,  (j  -  2,...,k)  except  that  i|^(l)  <.  fj(2)  for  every  j  -  2,...,k 

it  follows  from  the  above  result  for  n  -  2  that  (keeping  t|^(i)  -  i) 

n  k  .  n  k 

I  H  1.  I  H  SJ. 

i-1  j-1  *j(i)  i-1  j-1  *J( 


Hence,  the  original  permutations  could  not  be  optimal. 

Now  let  R(x^,...,x  )  be  any  real-valued  joint  cumulative 
probability  distribution  function.  For  our  purposes  we  want  to  prove 
the  following  extension. 


Theorem  1: 
n 


If  a|  (j  -  l,...,k),  then 

Ji  RUi . 1  R(aJ1(i)’  \(±) . \(i))’ 


where  il»^(i)  -  i  (i  -  l,...,n)  and  (j  -  2,...,k)  are  any 

permutations  of  1,  2,...,n. 


1.  This  extension  turns  out  to  be  a  rediscovery  of  a  special  case  of  a 
result  obtained  by  Lorentz  [2]. 


4 


a.<f  jPMW. 


Proof:  As  in  the  above  proof  we  need  only  prove  the  theorem  for 

n  ■  2.  Also,  without  loss  of  generality  we  can  then  take  a|  ■  0, 

a^  ■  1,  j  a  1 . k,  and  we  can  take  R  to  be  a  discrete 

distribution  function  with  mass  only  at  the  points  X^ . where 

(m  -  1,,..,2  )  are  all  the  points  consisting  of  the  k  coordinates 
which  are  0  or  1.  Let  Cm  denote  the  probability  mass  at  the 
point  X^.  For  any  ’f,2,*’**^k  we  can  wrde 


2  k 

"^(i) . \(1)>  -J,  Cm  \  «mj  <**  «))  ’ 


where 


gmJ  (a)  -  1,  1£  a  I  xmJ 


*>  0»  if  a  <  x 


where  x  .  is  the  j  coordinate  of  X  . 

mj  m 

This  is  so  since 


1  i 

n  8 -  °*  if  <  *,4  fol‘  at  least  one  * 


*»r  ^d) 


^j(i)  mj 


=  1»  ai  ^  > 


,  / j  \  >  x  .  for  all  j  =>  l,...,k. 
tj(i)  -  mj  J 


But  since  (a)  are  non-decreasing  functions  it  follows  from  the 
Hardy,  Littlewood,  and  Polya  result  for  products  that 


k  k  k  k 

n  g,(a^)+  n  gm,(ab>.  n  g.(ajj  ,,x)  +  n  g  ,  (a]  m) 

j»i  mJ  1  j»i  2  j.i  ^j(i>  j=1  mj  ipj  (2) 


K 

for  any  permutations  tj l>_,. . .  ,ik.  Then,  since  c  _>  0  (m  ■  1,...,2  ) 

z  cC  in 


(J  *  1 . k)  denote  the  set  of  numerical  values  of  the  k  components. 

1  k 

We  assume  that  R(a  ,  ...,a  )  Is  the  probability  that  the  system  will 

1  k 

perform  satisfactorily  (i.e.  R(a  , ...,a  )  is  the  reliability  of  the 

system)  where  R(a\...,ak)  as  a  function  of  a1,..., a*0  has  the 

properties  of  a  joint  cumulative  distribution  function.  For  example,  if 

the  system's  performance  depends  on  values  of  k  random  variables, 

Y^,...,Y^,  to  the  extent  that  the  system  performs  satisfactorily,  if 

and  only  if,  Y^  <_  ,  j  ■  1,2,.  ..,k,  then  if  Y^,...,Y^  have  joint 

lk  lk 

distribution  F(y^, . . .  ,y^) ,  then  R(a . a  )  -  F(a  , . . .  ,a  )  will 

have  the  assumed  interpretation  and  properties. 

Now  suppose  a|  <. .  .<_  a^  are  n  components  of  type  j 

(j  ■  1 . k) .  Then  n  systems  can  be  assembled  from  these  components. 

Let  N  denote  the  number  of  systems  that  perform  satisfactorily.  N 

is  a  random  variable  whose  distribution  will  depend  on  the  way  the  n 

systems  are  assembled.  A  direct  application  of  theorem  1  asserts: 

k-1 

Theorem  2:  Of  all  the  (n!)  different  ways  in  which  the  n  systems 
can  be  assembled,  EN  is  maximized  if  these  n  systems  have  reliabilities 
R(a^,. . .  ,a^)  (i  =  1 . n). 


i 

i  * 


6 


Proof.  For  any  assembly  defined  by  i|/^(i)  ■  i,  (1)  (j  -  2,...,n) 


2  i  2  k 

EN  -  )  R(a,  a, 

^1(1)  <l»2(i)*  *  i|»k(l) 


The  previous  result  holds. 

We  can  assert,  also,  the  following  result. 


1  k 

Theorem  3:  If  R(a^  (i)***'*®^  (i)^  —  1/2  for  every  1  and 

i^2> •  •  •  ><1^  then  the  same  assembly  that  maximizes  EN  minimizes  the 

variance  of  N. 


••  ^  k  1  k 

Proof.  Since  variance  N  "  I  R(aJ  (1), *  *  * t8J  <!)> (l-R(aJ(1) , .  •  •  »aJ  (i))) 

the  truth  of  the  theorem  need  only  be  established  for  n  *  2.  That  it 
is  true  for  n  ■  2  follows  from  the  following  lemma: 


Lemma:  If  1/2  £  P]_  i  £  P2  i  *■»  Pi  £  32  —  P2*  fnd  pl+p2  -  ql+q2 
then  p^l-p^  +  p2(l-p2)  <_  q^(l-q^)  +  q2(l-q2>. 

Proof:  On  subtracting  the  right  hand  member  from  the  left  we  get 

pi(1-pi>  +  p2u-p2)  "  qi^-qi)  ■  q2d-q2> 


2  2  2  2 
P1  +  p2  '  P1  "  p2  ‘  ql  "  q2  +  ql  +  q2 


ql  ”  P1  +  q2  ”  P2  +  ^P1  "  ql^  +  ^p2  “  q2^ 


L! 

i. 
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■  (qx  "  +  Px>  +  ^2  “  p2^q2  +  p2^  "  (ql  "  P1J  ”  (q2  ‘  P2J 

■  (qx  -  p^^  +  p-l  -  i]  +  (q2  -  P2)(q2  +  P2  ~  ^ 

<  (p2  ■  q2)(cii  +  px  -  D  +  U2  -  P2Hq2  +  p2  -  D 

■  (p2  “  q2)(ql  +  Pi  "  “  ^2  ‘  q2^  (q2  +  p2  - 

■  (p2  -  q2)(qi  +  Pi  "  ^2  ”  P2^ 

-  -(p2  -  q2)(q2  -  qx  +  p2  -  px>  1  °- 

Some  remarks  are  in  order: 

1.  Theorems  2  and  3  will  be  applicable  to  special  coherent 
structures,  namely,  those  where  the  reliability  function  is  not  only 
monotone  but  also  a  distribution  function. 

2.  It  should  be  stressed  that  the  exact  form  of  the  distribution 
function  need  not  be  known,  only  that  the  reliability  function  can  be 
assumed  to  be  in  the  form  of  some  distribution  function. 
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[3]  Another  Application  of  the  Hardy,  Littlewood  and  Polya  Result 
To  Reliability 

Consider  the  follov?ing  model:  A  stockpile  consists  of  n  units 
of  each  of  k  different  types  of  items.  Associated  with  each  item 

is  a  probability  that  the  item  will  perform  effectively.  We  denote 

i  fcH 

by  j  ■  l,...,k,  i  -  l,...,n,  the  probability  that  the  i 

unit  of  the  type  j  will  perform  effectively.  (This  item  will  be 

referred  to  as  the  item.)  These  probabilities  are  assumed  to  be 

independent  and  it  is  supposed  that 

P1  1  p2  1  1  pi*  J  “  1 . k. 

From  these  nk  items  we  must  construct  n  modules,  where  a  module 

consists  of  one  of  each  of  the  k  types  of  items.  We  say  that  a 

module  functions  if  all  of  its  k  elements  perform  effectively.  The 

problem  is  to  construct  the  modules  in  such  a  way  so  as  to  maximize 

the  probability  of  attaining  at  least  0  <_  r  <_  n  functioning  modules. 

Let  us  define  the  iC^  module  to  be  that  module  containing  the 

item  P^.  If  our  objective  was  to  maximize  the  expected  number  of 

functioning  modules,  then  the  Hardy,  Littlewood  and  Polya  result 

yields  the  solution.  Namely,  that  the  ith  module  consists  of  the 
12  k 

items  P^  P^,...,P^.  We  now  show  that  this  construction  not  only 
maximizes  the  expected  number  of  functioning  modules  but  it  also 
possesses  the  stronger  property  of  maximizing  the  probability  that 
the  number  of  functioning  modules  is  at  least  r,  for  each 
0  <  r  <  n. 
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Theorem  4:  The  probability  of  obtaining  at  least  r  functioning 
modules  is  maximized,  for  each  r  -  l,...,n,  by  letting  the  iC^ 
module  consist  of  the  items 


pi  pk  j  i 

a  *i  •  •  • 


Proof : 


Consider  any  arbitrary  construction  of  modules  -  call  it 


-  whose  1st  module  consists  of  the  items 


pi  p3 

ri»  r j  i  r j  >  •  • • *r j  • 

l2  1k 


where  1^  +  1,  and  suppose  that  the  sth  module  of  consists  of 

the  items 

pl  p^  P^ 

s’  Pl*  Pj3 . Pjk* 

Now  consider  a  different  construction  -  call  it  C ^  -  whose  1st  module 
consists  of  the  items 


pl  p^  p^  p*1 

1’  1’  .  £, 

3  k 


and  whose  sth  module  consists  of  the  items 


p1,  p,  ,  p3 . pk , 

s  V  n3  \ 


where 


ir  -  min(ir,jr),  my  ■  max(ir,jr),  r  -  3,...,k. 


Suppose  further  that  the  remaining  modules  at  are  identical  to 

those  at  C^.  The  probability  that  either  the  1st  or  the  sth  module 


of  functions  is 


p}.p?  -P l  . pJ  +  (i-p}-.p?  .P? . pJ  HPi.pJ.P?  . Pf  ) 

i  i2  13  H  1  M  K  8  1  U  Ji, 


2  3 


while  the  corresponding  probability  under  is 


p!.p?.p?  . P^  +  (i-p}.p?.p? . pJ  )(P*.P?  . P*  )  • 


1  1  1. 


'1  1  £. 


A  '  J  . *  / 

*k  8  i2  m3  \ 


As  the  number  of  the  other  n  -  2  modules  which  function  is 

stochastically  equal  under  or  Cj  it  follows  by  the  Hardy, 

Llttlewood  and  Polya  result  that  the  number  of  functioning  modules 

is  stochastically  larger  under  than  It  is  under  C^.  Hence  we 

need  only  consider  constructions  whose  first  module  contains  both  P^ 

2 

and  P^  Similarly  we  can  show  that  we  need  only  consider  constructions 

12  k 

whose  first  module  is  P^,  P^,...,P^.  Repeating  this  argument  on  the 
other  modules  completes  the  proof. 


We  make  the  following  remarks: 

1.  The  main  reason  we  have  considered  this  problem  is  that 
often  the  "total  system"  works  if,  and  only  if,  at  least  r  of  the 


modules  function. 

2.  If  we  suppose  that  k  ■  2  and  that  a  module  functions 

if  at  least  one  of  the  items  in  the  module  performs  effectively  then 

a  proof  similar  to  the  above  shows  that  the  construction  stochastically 

maximizing  the  number  of  functioning  components  is  the  one  whose  i^ 

1  2 

module  consists  of  the  items  P . ,  P  . . 

i’  n+l-i 
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